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FOREWORD 

The  area  of  Sonobuoy  performance  prediction  is  of  vital  importance  for  continued  im¬ 
provement  in  our  Antisubmarine  Warfare  (ASW)  capabilities.  As  part  of  the  Naval  Ocean¬ 
ographic  Office's  investment  in  ASW.  the  Environmental  Systems  Division  specializes  in 
on-scene  predictions.  This  report  identifies  the  major  on-scene  sonobuoy  performance 
prediction  programs  and  analyzes  the  characteristics  of  each.  This  information  should 
prove  a  valuable  contribution  to  future  fleet  improvements  in  on-scene  prediction. 


C.  H.  Bassett 
Captain,  USN 
Commanding  Officer 
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Ln  an  open  (.'in  iroimu'iU  wiiere  Uiere  is  little  or  no  inlonnatioii  aljout  tile  target,  the 
sinuilation  pi()>;i-anis  i;ive  \ai\in”  results.  The  data  in  this  report  show  that  these  pro;4riuns 
ei  aluaie  probatiiliti  of  deteelion  ({’(•ni  differently  for  the  same  set  of  patterns  and  input 
I  allies.  By  anali  zinii,  eumulatii  e  !’( if)  curves  there  are  general  relationships  that  can  lie  ob- 
sened  between  the  models:  (1)  TASDA  tends  to  be  conservative  in  comparison  with  AZ(  d 
and  SPAIM  at  low  figures  of  merit  since  il  assigns  a  POD  of  0  or  1  during  a  Monte  Carlo  tried. 
(2)  SPAM  is  tlie  most  optimistic  of  the  three  at  low  figaires  of  merit  since  it  uses  a  cumula¬ 
tive  time  window  luul  a  weighted  average  to  olitain  POD.  (3)  At  all  freciuencics  and  FOMs, 
AZOI  shows  small  variiuice  between  probabilit>'  values  for  each  of  the  patterns  since  tlie 
ladculation  of  POD  is  time  independent.  (4)  SPAM  is  more  conservative  in  its  POD  values 
at  high  FOMs  since  the  cumulative  time  window  builds  detection  probabilities.  (5)  Both 
A/.ol  and  TASDA  ha\e  a  sharp  lireakoff  point  for  FOM  below  which  POD  values  decline 
rapiiLLi  to  zero  liecause  ot  the  probability  of  detection  vs.  range  curves  used  in  each 
model.  Output  for  ADKPS  was  not  analyzed  since  the  number  of  patterns  it  evaluates  is 
limited  and  the  output  is  a  conditional  probability  of  detection.  Differences  between  com¬ 
puter  models  ncerl  to  lie  investigated  fiu'ther.  However,  further  investigation  should  be 
liasod  on  a  real  world  application  of  the  sonobuoy  programs. 


I. 


LNTHODLCTION 


Uecitlin^  on  a  "liest"  pattei’ii  is  not  lui  easy  task.  Ideally,  after  input  of  all  availatde 
data  and  e.xccution  of  a  sonobuoy  program,  program  output  provides  an  optimal  sonobuoy 
pattern.  There  are,  however,  manv  complex  decisions  to  make  in  specifying  a  pattern 
including: 


nhere  to  drop  each  sonobuoy, 

number  of  soncdjuoN  s  to  drop  in  :i  pattern, 

'(inoIiun\’  depth  settings, 


length  of  time  to  monitor  a  pattern,  and 
number  of  patterns  to  drop. 


-All  of  tlicse  decisions  ai-e  interrelated.  For  example,  tlie  numixir  of  patterns  depends  on 
tlie  size  of  the  search  area,  tlic  aircraft  on-station  time  and  the  availability'  of  other  air¬ 
craft.  .Additionally ,  tlie  niimljcr  of  Ixioy  s  required  for  a  particular  pattern  depends  on 
sonobuoy  depth  settings  and  monitoring  time  for  each  buoy.  Operational  constraints  and 
mission  goals  must  also  be  considered  to  provide  tlie  decision  malvcr  with  the  best  possible 
choices  for  each  decision. 

Four  of  tlie  tactical  sonobuoy  computer  programs  which  aid  in  these  decisions  were  in¬ 
vestigated  to  determine  tlieir  functional  effectiveness  as  an  aid  to  planning  search  patterns. 
The  programs  studied  were; 

1,  Tactical  Sonobuoy  Decision  Aid  (TASDA) 

2,  Algoritlim  for  Zone  Optimization  and  hwestigation  (AZOI) 

d.  Search  Pattern  Assessment  Model  (SPAM) 

-1.  Automated  Deployment  of  Sonobuoys  (ADEPS) 

rhis  report  proiides  a  description  of  the  simulation  techniques  used  in  each  model  in  Sec¬ 
tion  II.  Input  parameters  of  target  intelligence,  environmental  data,  buoy  data,  and  aircraft 
characteristics  are  compared  in  Section  III  followed  by  a  discussion  of  computer  output  from 
four  scenarios  in  Section  IT,  Sonoliuoy  pattern  definition  is  included  in  Appendix  A.  The 
propagation  loss  environment  is  described  in  Appendix  B.  Output  data  analyzed  for  this  re¬ 
port  are  included  in  Appendix  C’. 
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II. 


DFSCRIPTION  OF  COMPUTER  SIMULATION  MODELS 


The  tactical  programs  ADEPS,  TASDA,  AZOI  and  SPAM  provide  measures  of  effective¬ 
ness  for  sonobuoy  patterns.  ADEPS  uses  a  mathematical  technique  based  on  the  symmetric 
properties  of  certain  sonobuoy  fields  to  calculate  conditional  probability  of  detection.  The 
other  models,  TASDA,  AZOI  and  SPAM,  use  Monte  Carlo  simulation  techniques.  In  these 
latter  three  models  satistical  distributions  are  used  to  simulate  initial  target  location.  If 
target  position  and  speed  can  be  estimated,  a  normal  statistical  distribution  is  used.  Other¬ 
wise,  a  uniform  statistical  distribution  is  used  and  the  target  is  considered  equally  likely  to 
be  anwhere  in  a  search  area. 

Model  Description  for  ADEPS 

Two  search  methods  are  available  in  ADEPS  (Allison,  1970).  An  area  search  rou¬ 
tine  predicts  effectiveness  from  a  latticed  sonobuoy  pattern  with  equally  spaced  sonobuoys. 
This  routine  includes  the  following  possibilities  for  buoy  placement :  4-4,  2-2-2-2,  3-3-3, 
8-8,  5-6-5,  4-4-4-4,  6-6-6-6,  16-1(),  11-10-11,  8-8-8-8,  A  barrier  search  routine 
assumes  that  a  submarine  is  on  a  course  perpendicular  to  a  row  of  sonobuoys.  Only  the 
area  search  method  is  examined  in  this  report. 

In  the  ADEPS  urea  search  routine,  symmetric  characteristics  of  the  sonobuoy  field 
arc  used  to  divide  tlie  field  into  identical  subareas  as  shown  in  Figure  1.  A  lateral  range 
routine  computes  POD  versus  range  values  to  give  the  probability  that  a  sonobuoy  will 
detect  a  target  when  the  target  is  at  a  given  range  from  the  sonobuoy.  Ihdi\ldual  buoy 
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f  igure  1.  A  Symmetrical  Sul)section  of  a  4-4- i-1  Buoy  Field  in  ADFPS. 

Probability  of  Detection  is  calculated  for  a  s.\  inmetrical  sul:)- 
area  of  tlic  buoy  field  aiul  is  used  as  ii'dicator  of  effectiveness 
for  tile  entire  field. 

laU-r;d  i  ange  curves  are  combined  tu  determine  probability  of  detection  at  eqiudly  spaced 
points  ii.  a  subarea.  Conditional  probability-  (P(t))  is  cominited  from  the  sample  points  in  a 
.->ubarca  where  Uie  combinetl  proluibility  of  detection  is  greaUir  tluui  or  equal  to  0.5.  'llic 
avci  agc  distance  D  between  consecutive  detectioais  is  computed  by  siunpling  each  point 
sequcniialh’  fas  shown  in  Figure  5i  and  averaging  the  distances  between  IV)));!  0.5  sample 
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In  (<iimtion  1,  I’(ti  is  a  amditional  probability  of  dotcclion  in  a  suljarea  duriiA’  time  t. 

P(t>  -  1  -  o  rtt  0  (  i  ) 

Althouy,h  it  is  not  the  probability  that  a  suljmarine  will  be  detected  by  the  entire  sonobuov 
field,  it  is  a  measure  of  the  effectiveness  of  the  field.  The  final  product  of  the  model  is 
tlie  average  conditional  probability  of  detection  (P)  Riven  that  a  target  enters  a  symmetrical 
subsection  of  the  i)uoy  field,  ft  is  determined  by  averaging  the  values  for  J^(t)  over  the  on 
station  time  period. 

Model  Description  for  TASDA 

TASDA  operations  (NAIDC,  197())  are  divided  into  two  functional  objectives:  (1)  crea¬ 
tion  of  a  sonobuov  tactics  file,  and  (2)  application  of  Monte  Carlo  simulation  to  produce 
optimum  sonobuov  field  deployment  for  ASW  threats.  The  tactics  file  is  created  by  a 
geometry  definition  program  (GKOAIT),  TASDA  uses  Monte  Carlo  game  theory  to  stimu¬ 
late  100  random  tracks  of  a  target  through  a  sonobuoy  field  selected  by  GEOMT.  Aircraft 
time  on  station  is  divided  into  time  steps.  The  position  along  a  target  track  is  computed 
for  each  time  step.  If  signal  excess  is  greater  than  zero  for  two  consecutive  time  steps, 
a  d(‘tection  is  recordiKi.  The  frequency  of  detected  target  tracks  to  total  tracks  is  used 
as  an  estimate  of  the  expected  and  cumulative  probability  of  detection.  The  total  time 
a  target  is  detected  bv  one  or  more  sonobuoys  divided  by  the  number  of  unbroken  detec¬ 
tion  intervals  is  output  as  mean  holding  time.  The  time  of  first  detection  for  each  target 
track  is  recorded  and  the  total  of  these  times  divided  by  the  number  of  detected  tracks 
is  output  as  mean  time  to  first  detection. 


Moilel  l)esc'!'iptii)ti  for  AZ( 'I 

AZ(  il  (Hirnl):ium  niul  Dot-nv,  I!)""))  is  :i  Inc-lics  paclci^f’  wliich  inclLides  four  subprograms 
l.fKip,  (>P  ri,  WAPr,  and  (PAMI  .  i-xpected  pr()tjabili(\’ of  detecdion  d  l’nl))  b\-  a  sonoI)uoy 
field  is  computed  in  1,()()1>.  ()l>i'I  determines  the  pattern  spacing  which  pro\ifles  maximum 
f  lvip.  W'AI'I  calculations  incdude  cumulative  prol)al)ilit\- ol  detection.  f'.AM!  uses  lipf)!) 
obtained  1)\  !  i  >i  )p  to  estimate  ;m  optimum  buoy  spacinsi  for  a  yi\en  pattern  with  se.'ci’a! 
diflercnt  sonolniio  charaid  eristics,  .-^ludi  as  fip,ure  of  merit  ipt  I'M  i  or  bu()\'  depth. 


In  .\/.(  d,  rp()l)  is  estimated  bv  Ahnite  Cajdo  simulation.  A  lateral  rani;e  cur\  e  of 
prol)abilit\-  ol  detection  versus  signal  excess  is  given  bv  a  receiver  operating  curve  ITPtC'i 
\'.'hich  is  dcterniined  Irom  a  probabilii\  ol  lalse  idaian  (jf  pr  *  and  a  specific  s[)ecti'Lim 
aiialvvA'i'  set  lor  a  1  iw-miniitc'  intcgi'atioii  lime.  Signal  execs.-  '.yrsus  lauigc  is  detci'- 
mined  from  the  I'l 'M  and  the  propagation  loss  cur\'e.  The  two  tables  ol  \  ;ilues  are  eomliinf-d 
to  produce  P(d)  \  ersus  range  values.  Range  from  a  sonobur)\  is  deteririine<i  from  a  statis¬ 
tical  target  distribution,  neteclion  proluibilil v,  l’i,j»  1’^  sonobuot’  at  rruige  rj  from 

the  target  is  found  from  POD  versus  lauige  \;tlues  (I'igure  .'A.  In  equation  2,  is  the 

expelled  value  for  detection  when  the  target  is  at  a  parlieidar  position  /Xi.A'ji.  I'he  sul)- 
senpt  i  is  aji  index  used  to  number  target  positions  from  1  to  N. 


1  -  (1 


)  n  -  p. 


.(1 


(1  -  p. 


n  =  number  soroL'uovs 


PPOD  for  a  pattern  Is  eaieulatcd  bv  laliing  a  random  sample  of  N  laj.’get  ixjsitions.  Proba- 
bili!\  is  calculated  for  ettch  position  i  as  in  equation  2.  The  i:irot)abilities  are  summed  a,nd 
di\-ided  b\  N  to  calculate  KIHtI)  in  equation  A. 


Ll'ld': 


Z(x  ,Y  ) 


■ '  X 


,  f 


.Z(x  Y.,) 

H  i 


/n 


N  -  nuiiiher  of  target  posifionb 

(  'umulntivc  probaliMitv  of  detection  ((■|■>o^)|  is  also  ostimaterl  Ijy  Monte  (TtIo  simu¬ 
lation.  Tlie  matliematieal  lunction  used  in  tlie  computation  of  CPOD  is  lanomutl  cuul  gives 
Uie  probability'  ol  at  least  one  dcU-'ction  in  N  trials.  A  time  t^  is  chosen  randomly  from  .i 
liiiii.'  interv;il  I  using  a  i sindom  numf'er  generatoj'.  The  iiiterv.'il  1  is  defined  by  equation  A 
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I imger  be  on i t arc  i 


.\r  niH  ei'piint','  about  the  target  |x>si'jo,i  at  lime  ij  i.-  cah  ul  Ued  fnam  a  normal  statistical 
dpt  riljutioiv.  DcLoilion  probabililv  at  a  po-^ition  i.-  ■.  aJculaoPd  as  in  ejuation  2.  The  maxi¬ 
mum  expected  value  foi  det<'ction  bv  souo'aun  i  (figin'e  D  over  X  target  tracks  at  each 
time  tj  is  summcil  to  yield  eumutati'c  probaliiliti'  of  lielection  in  eijuation 
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The  CPOD  in  equation  5  is  an  estimate  for  the  lower  limit  of  cumulative  probability  of  de¬ 
tection  at  time  t.  A  large  number  of  target  tracks  must  be  generated  to  obtain  an  estimate 
of  CPOD  with  a  small  variance.  The  number  of  simulations  required  depends  on  the  signal 
excess  and  the  desired  accuracy  of  output  data.  The  probability  of  detecting  a  target  that 
has  not  yet  been  detected  is  calculated  from  CPOD  using  Bayes  theorem  of  probability. 


Figure  3.  A  Graph  of  POD  vs.  Bange  Values  Computed  in  the  AZOI  Tactics  Package. 


Figure  4.  Maximum  POD  Pj  j*  of  3  Target  Tracks  at  Time  tj  for  a  Sonobuoy  j 
in  AZOI. 
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Model  Description  i'or  SPAM 


The  ST’AM  program  computes  (’l’f)D  bv  M»)nte  Carlo  simulation  (Marin,  197(>).  The 
range  of  the  target  fi-om  a  sonobuo\'  is  found  l)\  a  stochastic  process.  Signal  excess  fSK) 
at  time  t  is  found  from  a  normal  distribution  with  mean  Si;(r>  and  a  standard  deviation 
of  f),  s,  10,  oi'  12  dB.  i'lie  probability  that  the  signal  excess  is  greater  than  zero  at  rringe 
r  is  gi\  en  by  equation  (i  (Marin,  l'J7()). 
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-  =  standard  deviation 
of  signal  excess 

To  use  equation  0  in  a  simulation  model,  it  is  neeessar\  to  introduce  a  time  element.  The 
expression  I'ittjt  is  introdueetl  as  tlie  proljabilitv  that  a  target  is  not  detected  at  r;uige  rj 
during  tiim  t;.  The  following  equations  are  assumed  for  the  time  element. 
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liquation  7a  states  lh;it  the  probalMIily  of  detection  I’jttj)  is  proportional  to  some  constant 
A.  plus  some  liighcr  order  terms.  The  higher  order  terms  in  eqitation  7u  can  l)e  neglocled 
for  a  small  time  period  t.  Equation  71>  states  tliat  the  target  is  either  detected  or  not 
detected  where  the  protjubility  Mj(0)  oi  no  detection  at  luiy  riuige  witli  no  time  elapsed  at 
tluit  r;ug>e  is  one.  Equation  7c  assumes  tliat  no  detection  at  time  t  and  no  detection  at  time 
t-  At  arc  independent  events.  The  assumptions  lead  to  an  e.\ponential  statistical  distrilxi- 
tion  for  I’jiti).  rhe  SPAM  pi'ogram  uses  the  eciualion  (Marin,  1977>) 

P.il,.!  -  1  -  e  for  t  0  '  fC 


in  a  simulation  technique  called  ihe  i‘umulalive  lime  window  (CTW).  The  lengtli  of  time  t^ 
ncecssar\  foi  a  sonolxio,\'  j  to  deteel  a  target  and  the  length  time  tj^  lietween  time  steps 
are  input.  Probabilities  I’jttjl  arc  computed  and  compared  with  random  numbers  to 
doU  iiiiine  wi'en  a  detection  occurs.  A  progrmn  tri;il  ends  eitlier  when  a  detection  occurs 
or  when  Uic  ;iircraft  on  station  time  cxpiics.  The  upper  limit  for  cunuilativc  probability  of 
dotoction  is  obtained  b\  comi lining  single  buo.y  probal>ilitics  independently  in  equation  1). 


•  1  -  -  I  .(t,;;  .  .  .  C  -  ,9) 

n  -=  number  of  sonobuoys 

A  luwr'i  limit  lor  probability  <^1  deteetioTi  at  lime  ij  is  olitaincd  assuming  complete 
drpendenc*'  as  in  cxiualion  pi. 
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n  =  number  ol  sonobuo\  s 

SPAM  outputs  cumulative  probatiility  slif^htly  less  tliim  hallway  lietween  the  two  extremes 
usins  equation  11. 


CPOD  --  .45  (CPOD  +  CPOD  ,  ) 
max  linn 
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Summary  of  Model  Descriptions 

I'ach  of  the  models  calculates  probability  of  detection  at  some  time  tj  in  order  to 
arrive  at  :ui  estimate  for  IIPOD.  .VDllT’S  calculates  the  average  conditional  probability 
in  a  symmetrical  subarea  of  a  buov  pattern.  TASDA  computes  proliability  of  detec¬ 
tion  by  one  or  more  buovs  (I’Dli  using  Monte  Carlo  simulation  of  100  target  tracks  and  a 
delect  (  Pdf)  -  1)  or  no  detect  (POD  -  Oi  analysis  for  ea^h  track.  uses  a  binomial 

distribution  function  to  calculate  l-  pol)  where  a  value  along  a  lime  track  is  chosen  randoni- 
1'. .  SP.AM  uses  Monte  Carlo  simulation  and  averages  T’jftj)  (previously  defined!  obtained  by 
assuming  dependence  of  sonoiniovs  with  Pjfti!  obtained  by  assuming  independence  of  sono- 
buovs. 

I  ach  of  the  models  estimates  cumulative  probabililv  of  detection.  .\DPPS  uses  an 
exivincntial  expression  to  estimate  CPOl).  'I'.A.SD.A  uses  a  detect-no  detert  Monte  Carlo 
mathematical  techni(|ue.  A/ol  uses  some  of  the  fiPOI)  calculations  to  obtain  an  estimate 
for  the  lower  limit  of  CPOI),  In  the  uniform  case  CI><)D  is  time  independent  and  is  equal 
to  PT’OD.  SPAAl  uses  a  cumulative  time  window  (CTW)  defined  by  Iq  and  t^  to  calculate 
CPCD. 


III. 


SC’KNAIUO  DKFlNri'K  )N 


Five  sonobuoy  iJiiUein.s  employed  b\  ASW  Uielieian.s  wore  ehoson  ko  daUi  ;uial\.si.s  in 
seenario  1;  the  eirelo  (cjllipse  in  tile  normal  ease),  1-1  1-1,  ehevi'on  (ehevi'on  skew  in 
tile  normal  ease),  5-()-o,  and  brnslitac'  patteiais.  I’heir  eonli};n rations  an;  ilelined  in 
AppeiulLx  A.  1'hc  same  inpul  tor  taiget  intc;llit;enei',  (nii  i  ronmental  data,  and  aireralt 
eharaeteristies  were  used  so  as  to  evaluate  the  elleetii etuess  ol  each  pattern.  ()ut|)ut 
w’as  generated  tor  one-half  hour,  one  hour,  and  tour  hours  ot  flight  time  in  all  eases 
preeeding  the  aircraft's  arrival  on-station.  'I’he  aiieraft  remained  on-station  four  liotirs 
in  all  cases  for  scenario  1.  The  circle,  ellipse,  1-4--1-4  and  ehe^r’oii  skew  patterns 
w’ere  evaluated  in  seenario  2.  I’he  on-stali<n)  time  was  extended  from  lour  to  1(1  hours 
in  the  second  scenario  and  output  wais  generated  lor  erne  hour  time  late.  In  seenario  .1, 
tlie  circle,  ellipse,  4-1-4-4  and  ehe\  ron  skewv  patterns  wi'i  e  searched  in  iiuarters  In 
foul’  planes  with  one  hour  time  iate  anil  four  hours  on  station.  Seenai’io  4  used  the 
same  target  data,  buoy  ilata,  and  aircraft  data  inpul  as  seenario  3,  however,  environ¬ 
mental  input  was  changed  for  otlier  comparison  purposes. 

Scenarios  1-3  were  based  on  tlie  following  environmental  input.  An  XBT  teinijcratiire 
profile  to  730  meters  recorded  0  July  1071  in  tlie  I'aeifie  Ocean  at  latitude  33“  3si'  North 
;uid  longitude  123“  30'  West  was  in(iut  into  tlie  Ditegrated  Command  ASW  Prediction  S\  stem 
(IC.\i>S)PBOFr.FN computer  prog-’".n.  This  program  merged  tlie  observed  deptli- 
temiieraturc  pairs  witli  a  deep  historical  profile  to  yield  a  total  temperature  profile  witli 
bottom  tleiitli  at  4000  meters.  The  program  tlien  calculated  a  sounti  speed  profile  from  tlie 
tot;J  temperature  profile  and  historicid  salinitv  values.  The  sound  speed  profile,  which 
e.xliibited  a  deptli  excess  indicative  of  convergence  zone  environment  (  I'igure  3),  was  input 
inU)  tlie  ICAPS  Fast  As.vniptotic  Coherent  riauismission  (F.ACT)  program  to  determine 
acoustic  propagation  loss.  By  varying  die  target  deptlis,  receiver  deptlis,  and  detection 
Irequcncics  inserted  into  the  FACT  program  nine  iiropagation  loss  (PL)  curves  were 
gcnerateil.  These  curves  appear  in  AppendLx  B. 

'I  he  environmenlal  input  for  scenario  4  was  a  Gulf  Stream  profile  from  latitude  39°  30' 
North  longitude  70°  29'  West  recorded  17  November  1909.  It  was  merged  with  1C.4P.S  deej) 
historv  data  to  proiluce  a  total  profile  with  bottom  depth  2400  meters  (Figure  Gi.  .Vn  .\t- 
I antic  slope  w.ater  profile  at  latitude  39°  33'  North,  longitude  71°  30'  West  recorded  18  Nov¬ 
ember  1909  was  mergeil  to  produce  a  total  profile  with  the  same  bottom  depth  (Figure  7i. 
Two  different  profiles  were  used  in  scenario  I  in  order  to  show  the  effect  of  van  ing  the 
environment  on  the  tactical  programs.  Since  scenaiio  4  uses  environmenlal  input  which 
is  different  from  the  other  '!  scenarios,  its  output  iloes  not  direcllv  compare  with  output 
from  the  first  three  scenarios. 


] 

i 


DEPTH  'JEL-XITY  sound  LtLXIT*'  PROFILE  LAT  3553N  LON  12559U 

OELXITY  M 'SEC  DATE  6  771 

0.0  1507.3  l'<60  1500  1530  15-40 


r- 


r 

Ui 

r> 

Oi 

r-, 


(7> 


00 


<I 

O 


i  E 

L 

yi  H 

^8 


G> 

'T 


a 

8 


®<S0®0©0<S'0;?'S>S>S0<S'0  0'S>SD’S> 

®^^O(VIC0C»'l/>'HCnC0'S)fXJ'S»'9O<S<5>'S''S)® 
r*'rr^ai)oi)0)®njTr-®»-<'0®®o®®®o 
^  ^  ^  ^  ru  I'n  T  'i)  OJ  "S  (XJ  l/)  ®  '»■ 

^*4  *-4  ^  nj  I'vj 


U) 


I’i^ine  ().  Atlantic  (lult  strenm  sounct  speed  pi'ofile  used  to  !::enernle  propajiat’on 
loss  curves  for  scenario  4,  fl'he  predominant  mode  of  sound  propa'f^a 
tion  i.s  In’  dnx'ct  path.) 


UELOCITY  SOUND  UELOCITV  PROFILE  LPT  3936N  L 

UELOCITV  n/SEC  DPTE  171169 


PHOC ESSING  R KQUm KRIKNTS 


r\ . 


All  tour  tuctic;il  sonobuoy  computer  progriuiis  require  buoy  input,  tart;et  input,  and 
eiulroiunenUil  input.  ADEPS  iloes  not  require  aircral't  input.  Some  input  values  are  preset 
in  each  pro^rmn  so  tliat  the  user  need  not  specify  a  long'  list  of  values.  Data  processing 
ir\puls  toi'  the  four  scenarios  considered  in  section  V  arc  given  in  Table  1. 

Buo\  Input 

P,uo\  input  wliich  define  sonobuot  patterns  for  evaluation  is  shown  in  Table  2.  In 
AI)E1*S  tlie  tottil  number  of  sonobuoys  input  must  be  a  multiple  of  eight.  Area  dimensions 
are  rc'ciuired  luid  arc  used  to  calculate  a  buoy  spacing.  The  buoys  are  equally  spaced  to 
co\ei  entire  area.  ADEPS  is  die  most  limited  program  in  terms  of  buoy  input,  Imt  the 
sMiimetric  patterns  allow  rapid  processing.  The  only  geometries  available  in  die  .ADEPS 
preset  pattern  packtige,  applicable  for  dais  study,  were  the  4-4-4-4  pattern  and  tlie  5-6-") 
pattern  witli  a  uniform  target  distriliution. 

r.ASDA  requires  the  number  of  sonobuoys,  Uioy  position  assignment,  the  number  of 
tlifferent  spacings,  and  buoy  spacing  limits  as  input.  Patterns  are  input  using  a  sonolxioy 
ixisition  plamier  chart.  Optionally,  die  user  may  input  buoy  positions  in  x-y  coordinates  to 
define  a  pattern.  GEOIilT  is  a  preprocessor  program  for  TASDA  and  creates  a  tactics  file 
of  up  to  20  sonobuoy  patterns  with  as  many  as  04  buoys  in  each  pattern.  This  preprocessing 
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h  .iiiifc  cotiiliinc'--  sio  cnsi'S  showni  in  Table  1  inlo  l(i2  runs  for  processing.  Kun  time  and 
Iniin  inpn!  t  i nu  are  ladLiCKl  since  1hio\  i^wsitions  are  simple  recalled  from  GKOMT  rather 
I  han  rout  i  nel\  eiitei-ed. 


Table  2.  Buoy  Input 


ADEPS 

TASDA 

AZOI 

SPAM 

Coordinates 

X 

X 

X 

Spacin.n 

X 

X 

X 

X 

Bu  o\  Trans  lo  r m  ation 

X 

X 

Preset  I'ailerns 

X 

tioometries 

X 

Aniliient  Noise  SDev 

X 

Di  AZOI  a  i)attern  {geometry  is  defined  by  specifying  a  set  of  x-y  coordinates  for  cacli 
l)uoy.  A  lactor  for  converting  relative  tlistance  to  nautical  miles  must  be  specified.  Tour 
t>  pes  of  transformations  may  change  the  basic  pattern:  x-spacing  expansion,  }  -spacing 
expansion,  pattern  shift  along  tlie  x-axis,  and  rotation  of  the  pattern  about  the  origin, 

Li  SPAM  the  user  specifies  tlie  number  of  sonobuoys  and  the  x-y  coordinates  in  nauti¬ 
cal  miles  for  each  ljuoj’.  Buoy  coordinates  are  input  in  the  order  in  which  tlic  buoys  are 
deploNcd.  As  the  buoys  ai’e  dropped,  tlie  user-specified  radio-frequency  (RF)  range  is 
usetl  to  determine  if  a  sonobuoy  can  be  monitored.  After  the  pattern  is  complete,  S1\\M 
simulates  monitoring  constraints.  The  user  specifies  the  average  number  of  buoys  that 
mi  aircraft  can  monitor  and  the  monitoring  time.  Using  the  value  for  the  average  number 
of  buoys,  the  program  randomly  selects  buoys  from  the  pattern  and  only  these  buo\  s  will 
be  used  to  determine  if  a  detection  occurs. 

Target  Input 

There  are  twenty  target  inputs  as  showm  in  Table  .3.  Li  ADEPS  only  four  of  tlie  inputs 
are  applicalile.  Target  speed  is  input.  The  area  search  option  corresponds  to  an  e\cnl\ 
spaced  position  input  where  tlie  target  remains  ui  tlie  search  area  at  fi.xed  positions,  llie 
liarrier  search  option  corresponds  to  tlie  univariate  normal  input  where  tlie  target  is 
assumed  to  progress  towards  a  iiai'rier  witli  a  target  track  perpendicular  to  tlie  harrier. 
Only  the  evenly  spaced  input  is  used  in  Uiis  study  to  seai’ch  an  area  with  .l-h-h  and  1-1-  1-1 
patterns  • 

Twelve  of  tlie  twenty  target  inputs  apply  to  TASDA.  Target  speed  is  required.  The 
target  type  must  be  specifieil  as  nuclear  or  conventional,  holding  or  trmisiting.  A  nuclear 
holding  target  was  used  for  this  study  (see  Table  1).  Target  location  may  lie  desci  ihed 
in  two  ways: 

1.  A  uniforni  time-dependent  densify  function  assumes  that  the  target  deiisitx  is 
random  over  a  designated  rectangle  (Figaire  S). 


1.3 


2.  A  Bessel  normal  distribution  function  is  j^eneruted  from  an  initial  tjivariate 
normal  density  function  that  is  expanded  at  a  rate  equal  to  input  target  speed 
(Figure  9). 

Hotli  target  location  functions  were  used  in  tliis  study. 


Table  3.  Target  biput 


Target  Movement 

ADD  PS 

TASDA 

“  '  ""  1 

AZOI 

SPAM 

Course 

X 

X 

Course  Standard  Deviation 

X 

Course  Limits 

X 

Speed 

X 

X 

X 

X 

Speed  standard  Deviation 

X 

X 

t’clocitv  \’ector 

X 

I  Conventional 

t 

X 

.\uelear 

X 

Holding 

X 

Transiting 

X 

Distiuice  Traveled 

X 

Snorkel  Cycle 

X 

t  arget  l  ocation  Standard  Deviation 

X 

X 

Initial  Target  l.ocation  DeiisiD  I'unctions 

iivi.'iilv  spaeetl  positions 

X 

Ciii  forni  Time-tlependent 

X 

1  iiitoiin  Tinie-indeiiendenl 

X 

X 

ltd'  ariale  nornial 

X 

X  ! 

'  1  i- \  iionvKil 

1 

X 

X 

i;i\  ariate  normal 

X 

X 

J  iine-di:pi  iiilent  bivariate  normal 

X 

i 

\/u[  uses  sexen  ol  iwent',  Laiget  iitpuls.  Tiu’get  speed  and  speed  standard  deviation 
are  i:  pni  (lire(’tlv  '>r  in  teriii.'  el  ,i  .-ipei'<l  eovai'iaiiee  matrix.  I  tie  following  four  target  loea- 
lien  :  :i  funetiDiis  an  used  to  |iredict  initial  target  liX'ation  for  Monte  Carlo  simulation: 


RANDOM  TARGET  TRACK 


I'ip,T.ire  S.  Cniformh  nistiibuted  Initial  Target  Area  in  TASDA.  The 

target  location  is  random  over  a  designated  area.  The  same 
type  of  target  area  definition  is  available  in  AZOI  and  SPAM. 
However,  the  statistical  distribution  is  time-dependent  for 
lASDA. 


Figure  0.  Normally  Distributed  Initial  Target  Area  in  TASDA.  The  ini¬ 
tial  target  location  is  normally  distributed  around  (XPT,  YPT) 
with  standard  deviation  SGMAX  in  the  x-direction  and  SGMAY 
in  the  y  direction.  The  same  type  of  area  definition  is  avail¬ 
able  in  AZOI  and  SPAM.  However,  the  statistical  distribution 
is  bivariate  normal  for  SPAM. 


1.  A  unit’onn  tim<‘-indoiK;ndcnt  densits  function  assumes  that  tarf^et  clensit>'  is  con- 
st;mt  over  a  designaUid  area. 

2.  A  Bessel  normal  density  function  is  generated  from  an  initial  bivariate  normal 
tlensiU  function  ;uid  exi);unls  at  a  rate  equal  to  tile  best  estimate  of  tai'get  veloeitA 

A  time -dependent  bivariate  normal  density  function  predicts  a  target  position 
when  target  Iracl-dng  lias  been  suceessful.  This  function  describes  a  target  being 
trackeil  when  if  is  imperative  to  re-establish  target  contact, 

1.  A  bivai'iate  normal  distriliution  function  assumes  tliat  a  target  is  normally  dis¬ 
tributed  in  mi  area  generated  from  a  single  DIFAK  buoy  fix,  a  SOSl'S  conuict,  or 
radar  contact.  This  distribution  does  not  ex])and  with  time. 

rhe  fii’.sL  mid  second  funetions  were  used  for  uniforni  mid  normal  cases  in  this  study 
(Table  1). 

.SPAM  uses  seven  of  the  twent\  inputs  listed  in  Table  3.  Bequired  input  is  speed, 
speed  stmulard  deviation,  course  mul  course  stmidard  deviation.  Target  location  cmi  be 
assigned  using  tliree  distribution  functions: 

1.  A  uniform  time -independent  distribution  function  witli  the  propert}'  that  a  target 
is  equally  likely  to  be  at  any  point  within  a  bearing  box  of  length  L  and  half¬ 
width  H. 

2.  A  bivariate  normal  distribution  function  which  has  a  user  specified  mean  mid 
standard  deviation. 

:i.  A  univariate  normal  distribution  function  which  has  a  normal  distribution  in  an 
\  direction  mid  a  uniform  distribution  ui  a  y  direction. 

I'unctiun.-,  1  mid  2  were  used  for  data  mialysis  in  tliis  stud\’. 

/va  usti'j  Liput 

.'.rnnstie  input  for  e.'icli  program  is  listed  in  Tatde  -1.  The  acoustic  inpu'. ,  source 
level  (SI,),  mnbient  noise  (AX),  mid  recognition  dilierential  (RD)  are  usee  to  liotermine 

IS  i;.  tifjiiation  12. 

ton  M  Rp  ;:2 


Prob.ninlit>  of  detcetioi:  is  del'inod  ti,«  be  >).7,  when  signal  excess  (SE)  in  equation  13  is  zero 
tor  pi (ipag  xtion  less  (PI.)  at  '’ange  R. 

r  ^  FPd  -  ni  wO  1 1: 


p- 


Lii  |uii:,;r;un  a  iiioliabiliU  CLuve  is  used  to  assign  probability  values  to  values  of  signal 
exc  i'ss.  \M\eii  sigiuil  excess  is  greater  tluui  zero,  tlie  probability  values  will  be  greater 
Ulan  ii.  j.  \Mu'ii  signal  cnxeess  is  less  than  zero  proba))ility  values  will  t:)e  less  than  o.  b. 


Tal)le  1.  Acoustic  Input 


AD  EPS 

TASDA 

AZOI 

SPAM 

I're.p  I  (,ss  (PI.) 

X 

X 

X 

X 

S'  ■ina  (  !  .c  VI  ■] 

X 

\nd)ienL  Nc.iise 

X 

l!ccognili(.>n  IlitlerenLi;J 

X 

J'igui'c  of  .Merit 

X 

X 

X 

.Multiple  FO.M 

X 

X 

X 

lit  If  ('ur\  (.1 

X 

Multiitle  I’L 

X 

X 

X 

C on V e r g e  1 1  e e  Zone s 

X 

X 

AD)  I’S  considers  propagation  loss  (PI.)  curves,  source  level  (SL),  ambient  noise  (AN), 
and  recognition  differential  (it!)),  but  ADKPS  does  not  have  a  pro\ision  for  time  late.  As 
a  result  ADl  l’S  could  onh  l)e  applied  for  34  of  the  810  possible  cases.  There  is  no  pro- 
■rision  for  multiple  fre((uencic.s  or  multiple  FOIMs  and  a  separate  run  was  required  for  each 
case.  In  order  to  reduce  tlie  number  of  runs,  an  FOM  of  75  was  selected  with  a  frequency 
nf  :;nfi  ilz,  t  ighteen  comiuiter  runs  were  processed  for  scenario  1. 

The  T.ASDA  progratn  can  accept  up  to  tliree  propagation  loss  curves  and  up  to  four 
FOAIs.  These  tw'o  features  reduced  Uie  number  of  runs  required  from  102  to  18.  The  user 
may  input  convergence  zone  ixmges  and  widths  instead  of  propagation  loss  curves. 

In  AZCH,  a  receiver  operating  characteristic  (HOC)  curve  is  required  input.  I'his 
cun  e  gi\  es  pi'ol)rilnlit\  values  ns  a  function  of  signal  excess.  Other  inputs  are  a  table  of 
propagation  loss  versus  range  and  a  figure  of  merit.  Tw'o  subprograms  of  .\Z01  were 
modified  in  Ma\  l!)7s  to  process  multiple  FOMs,  multiple  propagation  loss  curves,  and 
multiple  patterns,  fhe  number  of  required  data  sets  for  input  w’ere  therebv  reduced  from 
sff)  to  eight. 


''I’Wl  ;k\' t'pt  s  ( UK  Jin  )paL',:ii  ii>;i  In.-  Ii!i  ■'.ii!'  I  r( '  i.kmic',  -iiiiI  i  cceix  e  i'-- --iiii  rc  t  ilM  1^- 
si;i  >  il(pih  i(k‘iU  i  ru':ii  iMii  luimin  i-..  .\ii\  ikhuIk  r  ol  pi'dpau.iiioii  |.iss  KiiiKt'-  iti.'r,  iiu  infiiM 

In  liu  lili  vith  :mi\  iiunilu  r  cil  I  r<.'<|iu  lu  ii  i'l  juloniat  icall  \  ffhi''-  iliionuli 

lh<-  iiipul  propapal  ion  lo-.a  tilt-  i..;  i1k  .app'i'i.i'ia.ili  |l|■n|)a^,•^^ioll  losa  cuikc.  I  la  main  pi’o- 
piaiv.  '>I’AM  w  a>  nioditinl  in  p!nml)L-i-  laTT  I'or  nin  ll  ii)i<.'  I  (  >M ,  imiltipU'  lime  late,  and 
m'.iUipK  propaKalion  loss  input.  I  lu'  nunil'ef  of  oxeeidioiis  was  reduced  Ifom  s]()  p,  ci"ht. 

\  I n  1 'pap.U ioi I  Ins-  eui'Ni'  is  etas--iiui!  a.s  -liallo". -shal lo','.  (.S-Si  lor  :i  shallo'.'.  iari-:,el  and 
I,  n(i\ti,  I  i<  1  p-dt,  t  p  I  I )' I  )i  tor  a  deep  '  'rpel  .and  ra'cei'.ei',  oi'  sh.al  iou-deep  i  S- ! )  I  I'oreitliei 
,1  -la'lo,,  tar;ael  or  rei'ci  vt' r  or,  output. 

.\iia  ladl  input 

.XDirL'S  rloc-S  not  eumsider  :Lireiaili  paiaunuLei  s  cxeo[iL  ior  monitoi'  '  lepi  adalion  a.s  .shown 
in  i  ;Ll.)le  ).  i'ur  TASDA,  Uie  Lli:j,hl iJuUi  is  dcsci  iljcd  oial\  Lu  Uie  exlcnL  (ji  sijecih  in^  Uie 
■  iap-ed  Lime  licioro  ;u'ri\aLl  on  slalion.  I'lie  r.alio  liequenc\  (liF)  1';UV4C  is  Lite  nnuximum 
.illowable  di.spmce  IjcLwcer,  Llu'  aireraiL  ;uid  a  sonohuox  wliicl)  pcrniiL.s  I'atlio  rccr'pLion. 
i  or  r.'\Sl).\.  th('  PusiLion  ol  tJie  .lirerail  is  ii.xed  al  Uie  eenLcr  ot  tJie  iielil  ior  Uie  sinuilaLion. 
\/i  i|  uses  time  late,  airerat't  ■\adocil\,  and  on-stati(,>n  time  to  ealcuhtte  fl't'l).  .SI’.AM  in- 
elucies  lim,.  late  :ind  tinie  on  staiioii  as  •,\ell  ,a.s  aircr.ali  Ki'ound  speed.  I'he  iis£  i'  spiecified 
-P’oui’.il  i-pi.a'd  in  SP.XM  deteniiines  the  Lime  at  whicli  each  buo.\  comes  up.  f)nl'.  Uiose  buo'.  s 
.iLirii'-  a  usv'r-.spi'eiiied  di'  ranee  tire  inoniloi'c,-d.  .\iter  all  l>uo,\'.s  are  in  Die  wtiLer,  Lite  user 
mil!  a  sonoi)UO\  .noniU.u'in^  selieme'  Liial  r-jllects  i!  F  e(;nsLi’ainLs. 


rablc  Ail  cral't  Liitiil 
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'  ;\I>1.1>S  j  i'ASDA 
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^  Hi'  Hanpe 
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1  '.In  litoiiiie,  Sequence 

I 

1  X 

( 'oivl  .S'leral 

'lonil or  l>,'K|-nd;ni,,ii 

\ 

1 

1  X 

_ L _ 1 _ 

i 

sn  i  :  m  :  ol  ;\|oi  lel  input 

i.a-  ■  oi  'Ja  ;  .1  Op  i uum-.  use.-  diilei  eiii  inaUiemalietu  I'ai  inulaLions  Lo  rleseribc  a  laelieal 
-it'  alio  ll.  'ipeii  (.aauu:  'au  iroiijnei'L.  Tlic  ionniil.iUon:-  are  based  on  as.su  inptions  alxul 
Uie  arianea  '  1  Uie  I  cnM,  Lin  ajnourit,  ol  .signal  e.\ce.s.s  I'eeos.sary  to  detect  a  lai  e<;L,  ;uid  U'.e 
'ependc  CO  (('1  inde[wndencoi  ol  .sonobu.jv  eolloetivcu  contributiae;  to  a  detection. 


1" 


'.i>!  '  11-  '1  pnKi'itis.  I  iic  presc-f  p;ickaf;e  caiinof  l)e  modified.  TASJ)A  also 

('  paiU  in  package,  Inii  ihe  package  can  l)c  modified.  AZoi  and  SPAM  have  no 

piH  sei  pattei-n  files.  ADI  1>S  has  a  normal  initial  target  location  distribution.  TASDA, 

\'/‘  'I  and  SPAM  all  have  oi'tions  for  uniform  target  location  distributions  and  normal  loca- 
'  ini'!  dist  ri  but  ions. 

\P  miidels  input  a  panpagation  loss  curve.  TASDA  may  input  propagation  loss  curves 
at  1'"  ,!,,  different  deptlis.  Therefoie,  one  sonobuoy  pattern  may  have  sonobuoys  placed  at 
dii:.  ri'iii  depths.  A'/(  d  inputs  a  receiver  operating  characteristic  tltoC)  curve  which  assigns 
'  I 'lobai'ilily  to  \  allies  of  SI  .  In  ADKPS  and  SPAM  this  curve  has  a  normal  density  distri- 
a.tion  witli  proliaihlitv  e(|  ;al  to  0.5  at  the  mean  of  the  normal  distribution.  In  TASDA 
ctin  tipuive  prnbability  i-  determined  from  a  Pician  density  Distribution.  'i'ASDA  assumes 
all  iioi)uo'\  s  are  :tcti\  atod  :it  tlie  same  time.  In  SPAM,  the  user  specified  ground  speed 
d»nei  mines  tiie  time  each  buov  comes  up. 

ANAl.YSIS  (  il'  Pf  SIT.'l'S 

I  wo  t\pes  of  program  outnu'  are  compared  by  this  study:  expected  probability  of  de- 
1  ection  i  i  pt  )Di  and  cumulative  ,>robability  of  detection  (CPOD).  KPOD  follows  the  statistical 
definitions  for  expectation.  CPdD  is  computed  from  a  probability  density  function  which 
follows  the  statistical  definition  of  a  density  function.  Each  program  uses  different  assump¬ 
tions  to  arrive  at  i:P(iD  and  different  density  functions  to  calculate  CPOD,  Since  methods 
of  calculations  are  different  in  each  program  the  output  have  different  numerical  values. 

The  anah  sis  in  this  seelion  examines  mathematical  techniques  to  explain  the  CPOD  and 
I'  l’i  )D  numerical  output. 

A  list  of  all  program  output  is  given  in  Table  6.  Although  CPOD  and  EPOD  para¬ 
meters  in  Table  5  are  listed  for  TASDA,  AZOI,  and  SPAM,  numerical  output  for  these 
measures  of  pattern  effectiveness  are  not  the  same.  The  output  values  were  analyzed  in 
reference  to  four  scenarios.  Three  target-receiver  depth  combinations,  shallow-shallow 
(S-Sl,  shallow-deep  (S-D),  and  deep-deep  (D-D),  were  used  to  generate  propagation  loss 
cuiwes  for  the  fir.st  3  scenarios.  Probability  of  detection  was  predicted  for  three  FOMs  and 
for  a  low  frequency  range,  middle  frequency  range,  and  a  high  frequency  range  in  scenario 
1,  2  and  .3.  Scenario  I  uses  two  target-receiver  depth  combinations,  two  frequencies  and 
one  PfiM, 

Scenario  1 

Expected  probability  of  detection  was  calculated  for  five  patterns  with  a  uniform  and 
normal  initial  target  location  distribution.  Time  late  was  not  a  factor  in  the  uniform  case 
because  that  parameter  did  not  change  the  initial  target  location  distribution  for  AZOI  and 
SPAM.  Therefore,  values  for  1  and  4  hrs  time  late  have  been  omitted  and  only  values  for 
30  min  time  late  appear  in  Appendix  C.  Cumulative  probability  of  detection  was  analyzed 
using  the  method  shown  in  Figure  C-1.  Esing  this  method  cumulative  probability  of  detec¬ 
tion  was  converted  to  a  single  eventual  probability  of  detection  value  for  comparison  pur¬ 
poses.  Both  the  eventual  POD  values  and  the  original  CPOD  curves  appear  in  AppendLx  C 
for  TASDA,  AZOI,  and  SPAM.  The  18  CPOD  curves  that  were  generated  for  ADEPS  are 
shown  in  Figure  10  with  the  corresponding  curves  for  TASDA,  AZOI  and  SPAM.  Since 
ADt;ps  could  process  only  two  of  the  patterns  the  results  were  not  analyzed  for  inclusion 
in  Table  7.  However,  the  relalionsliip  between  ADEPS  output  and  output  from  the  other 

ly 


tln'ff  models  c;m  lie  detluee<l  from  I'i.uiii'e  10.  (me  characteristic  of  the  Allld'S  curves  is 
that  ihev  al\va\  s  rise  to  100'  I’til). 


Table  G.  Measures  ol  Kffoctivcncss 
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Pnilorm  (  ase 


cimiulutivc  probability’  as  a  jfu  action  of  time  on  station 
expected  probability  of  detection  on  one  or  more  sonobuoys 
probability  of  detection  on  two  or  more  sonobuoys 
l)robabilit>'  of  detection  on  tliree  or  more  sonobuoys 
mean  holding  time  on  one  or  more  sonobuoys 
me;in  holding  time  on  two  or  more  sonobuoys 
me:ui  I'.ulding  time  on  three  or  more  sonobuoys 
memi  Lime  to  first  detection 

avei  a.ge  conditional  probabilitv  in  a  section  where  tlie  submarine 
is  assumed  to  be  a  random  target 

eumulative  probiUjility  tliat  a  submarine  has  entered  part  of  a 
liuu>  field  where  signal  excess  SE  is  greater  tlian  /.ero 
stmidard  deviation  of  POD 

probability  tliat  a  target  will  be  dfdected  given  that  it  has  not 
hecji  detected  in  tlie  past 


[;P(  1)  values  for  the  uniform  case  were  averaged  over  three  FOMs  and  three  frequen¬ 
cies  and  ihe  averages  are  present  in  Table  7.  For  a  shallow  target  and  receiver,  TASDA 
diuvvs  'l'.u  1-1-  I-  I  tvatlciTi  as  most  I'ffective  over  the  frequencies  ;uid  FOMs  with  a  value 
if  'll.;'  !  p(>D,  l  or  the  shallow-deep  target -receiver  combination,  TASDA  shows  the 

chevron  ixittem  as  most  effedi'  e.  For  a  deep-deep  target  receiver  propagation  loss 
curve,  CASDA  selected  the  clnn  ron  pattern,  A7.(  n  selected  the  circle  pattern  and  SP.MM 


1  I 


•  lASDA 
AAZOt 
OS  I ’AM 

5-  6-5  PATTERN 


Figure  10.  I’robability  of  Detection  \ersiis  Flapsed  I’ime  on  Station  for 
Scenario  1  with  a  Fniform  Target  Location  Dislril)ution. 
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selected  the  4-4-4-4  pattern  as  most  effective.  Note  that  there  is  very  little  agreement  in 
pattern  selection.  This  is  because  of  the  different  modeling  in  each  program.  Since  the 
patterns  were  defined  at  a  near  optimum  sonobuoy  spacing,  the  models  give  about  the  same 
probability  averages  for  each  pattern  they  evaluate.  However,  the  difference  in  average 
probability  of  detection  from  one  model  to  another  is  very  noticeable.  Differences  in  model 
output  is  as  much  as  34*^7  predicting  the  performance  of  the  same  pattern.  Since  there 
was  not  much  difference  in  i:pOD  between  the  five  patterns  using  one  model,  the  number  of 
patterns  consideretl  was  reduced  to  two.  Only  the  circle  and  4-4-4-4  pattern  were  pro¬ 
cessed  in  subseciuent  scenarios. 
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Nni'inal  I'ase 


l  ime  late  is  a  factor  for  1-.P()1)  in  the  normal  case.  The  time  it  takes  an  aircraft  to 
arrive  on  station  changes  the  distiibution  for  initial  target  location.  At  time  late  equal 
zero  I'ASDA,  AZoi,  and  sl’A.M  all  have  the  same  normtil  target  location  distribution  as 
shown  in  Tigiire  11.  At  .‘10  ir.in  lime  late,  the  norm.al  distribution  spreads  out  as  shoum 
in  I'igure  111  and  Mgcire  1,1.  .After  .‘10  min  the  probability  that  the  target  is  at  the  original 
estim;Ued  position  is  smaller.  Since  tlie  target  location  distribution  does  change,  the 
eventual  probabilitv  of  detection  v.alues  are  different  for  each  lime  late.  These  values  are 
given  in  .Appendix  II  under  Noimttl  Case,  in  this  study  normal  case  output  Cl’OD  values 
were  av  er;ige<l  over  tlii’ee  T(>Als  and  three  frevpiencies  and  over  each  time  late.  Table  7 
>h(A\s  .vpi'uged  values.  1  ho  three  sets  of  output  show  that  .AZOI  is  consistentlj'  the  most 
optimistic  of  the  three  programs.  This  is  a  change  from  the  uniform  case  where  SPAAl 
was  the  most  o|vtimi.stic.  .As  in  the  uniform  case,  the  averaged  values  were  very  close 
foi-  difierent  patterns  using  the  same  model.  The  ellipse  and  chevron  skew  patterns  were 
Ijrocessevi  suiiseviuent  seenaiios. 


.scenario  H 

(■'umdativ  e  pisjbabilitv  of  detection  was  calculated  for  tvv’o  patterns  with  a  uniform  and 
normal  initial  target  location  distribution.  I'he  first  llight  arrives  on  station  after  one 
hour  lime  late  ;uh1  remains  on  st.ation  for  four  hours.  A  standard  deviation  in  target  speed 
i^  introductHl.  Th.e  last  th.ree  llights  arrive  on  station  with  no  time  late  and  remain  on  sta¬ 
tion  foi'  tour  hours,  this  scenario  is  .simulated  bv  a  comitnter  ntn  with  one  hour  time  late 
.'Jill  pi  liou rs  on  station. 

Tile  values  in  i  able  s  were  obt.unecl  by  examining  Cr>ol)  curves  and  attaching  a  numeri- 
t  .al  V  ..lue  to  e:ich  cui  v  e.  In  order  to  arrive  at  a  numerical  value  the  slope  of  the  curve  and 
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Percent  TroPabilitv  of  Detection  for  .Scenario  2 
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FiSini  e  12.  Normal  Target  Location  Distribution  for  I’ASDA 
and  AZOI  at  30  min  Time  Late. 
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Figure  13.  Normal  Target  Location  Distiibution  for  SPAM 
at  30  min  Time  Late. 


its  maximiini  \  aku*  sliouki  he  considered.  A  lunnerieal  \:ilue  linscd  on  nren  under  n  ciirxc 
used  since  area  under  a  cur\e  depends  on  rise  time  and  cur\e  heiahl.  IA>r  exampk,  in 
l'i'-;T.ii  e  11,  C'ina  es  A  and  15  rise  h;  (he  same  pi ol)al)ililv  of  detection  etjual  to  L'O.s,  ll(i',ve\er, 
cur\e  A  rises  to  that  I’fil)  faster  than  ciir\e  15  and  has  a  hi^;!ier  numerical  indicatoi'  at'aciied. 
In  I'iaaire  In,  cur\e  ('  rises  as  fast  as  cur\-e  A,  Init  curve  ('  rises  to  a  -^mallf  r  I’OI)  maxi¬ 
mum.  Therefore,  the  lactic  nssncialed  with  curve  \  has  a  higher  measura-  of  elteci  iveiu-T^ 
than  curve  ('  and  thus,  has  a  Itipher  numerif-al  indicator. 

The  area  under  the  t'Pt)!)  curxes  w;>s  estimated  h\  \isuall\'  inspecting  the  curves.  The 
method  used  to  calculate  the  CIHil)  curve  indicator  '.allies  is  slujwn  in  I'iKiire  (.  -2  in  ,\piten- 
di\  ('. 

I  nil'orm  Case 

Cl’i  'O  indicator  \'alue.s  for  llu'  uniform  case  were  at  eraftetl  over  three  ft  isls  and  three 
tre<(uenci(  s,  and  the  ax'erttses  are  pia.sentetl  i'l  Tahl<‘  s.  Tlie  indr  sitor  \  allies  lor  e.ich 
l  u  r\ f  are  shown  in  Appendi.x  C.  In  this  seenario  onh'  two  patterns  were  considered  ;ind 
computer  out|)ut  showeil  onlv  partial  agreement.  However,  it  is  important  to  remem.her 
that  in  the  uniform  ctise  initial  target  [Xisition  i.s  nnknown,  I'or  all  three  models  indicator 
i.alues  fell  with.in  ti  few  percentage  values  of  eticli  other  for  the  1-1-4- 1  and  circle  paiierns. 
However,  the  difference  lieUveen  motlels  is  ;is  much  as  ;>7  in  the  14-1)  case  and  i,-  poiicc- 
ahle  in  all  ettses. 

Normal  Case 

Cl’i  lO  indictitor  tailues  for  ih.e  normal  crisc  were  ;r- or  aged  ot'cr  three  TCMs  and  tfiree 
fre(|uencic's  and  the  averages  are  siiown  in  i'alile  Thtj  indicator  \  trlucs  are  sTfOwn  in 
.Apiicndi  -  (  .  In  this  case,  ini'ial  target  location  ijoin's  are  clustered  towartl  the  center  of 
the  opi*:  aiion  area  for  I  .\s|)  \,  h/td  ,ind  Sih\.M,  (hne  would  expect  these  three  models  to 
s.'ahuiu  Ct-  nai'etiis  i  mi  Icrh  since  there  i-'  now  .some  control  on  target  location.  I  ;ihle  8 

■  iaw'  pi-()li  ibi !it\  ot  deteciioa  rt-main  not iccaltiv  different.  Note  that  '.he 

:  ,\om  a,  ■ ' i ! ;  oreri i-  ii  'iu  S- .S  iepti]  w  itii  the  ellipse  whereas  in  the  tiiiiforni 

'  -  :  !■  a  '  I’l  la.  :'l)  il  ei  er.i  <  '.''as  _’-l  M-i  11, e  eircle  am  i.lic  S -S  depth,  i'his  relatio)  ...-lup 

■  '  ■h'fcrem'i-s  can  be  obs..r\ed  for  ho'h  paitcrns  at  all  depths.  The  long  on  station 

C'  ■  '  b  hoc  I  s  api  i  a  I’s  lo  ini .  c  an  aflecl  on  t  lie 'urge  differences. 


'  "  ‘  c  i  p  1  ob.ii  .il  il  \  -)i  i!eiceiii,ii  inMsuaJ  Was  calculaUil  for  t'.'o  paitcrns  tvitl;  a  uni- 

1  I,  m.'i]  iCi.al  ;  'fget  I  oc  il  ion  '  ii  si  r’i  ’  r  u  ior .  T.itir  llCghis  arrive  on  station  with  I 

'  ■  «.  ai  d  la.ni.'uii  ua  station  for  torn'  houi-s.  !  ach  High'  searches  a  (juadrant  of 

'  ■  '  b(  )d.  Siam.  Mil  siiudriov  pattern.-  ar'  .'  mrueti-ic  with  rcspiect  to  the  center 

■1  •'  I  '  .  ,  1  acii  .|ii,iilraji'.  has  an  i'h.niical  so.iob'jo\  p,e<vn..4rv.  Therefore,  only  computer 
-iu.  I  1.0  on  o'm  iiiidrani  were  necessar'..  'r.sie.ati  of  obtaining  C  I'Clt  indicator  t'alues 
'  ;  .  .  '!>'  '  linn  (.nr  es.  a,,  ra',e  '.alues  were  oblaine.il  i)\  diiiding  the  Cl’'  >1)  curtes 

'o'lini  .1(0-  '  ,  ,  i,  i'i  1,1,  .  . '  of  Apiicndi.x  '  .  i '  I't  d  )  in'licator  values  tor  tlie  iini- 

!o;n.  ■:  o  nr'  I  ease  wn  re  j-  |._ed  o'rcr  three  C'  Al.o  and  three  fretiucncies  and  are 
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Fifixire  M.  ^  jrves  A  and  B:  Numerical  Value  for 
Iritfem  Selection  by  Rise  Time. 
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Fisiire  l.l.  Curves  A  and  C:  Numei’icjil  Value  for 
pattern  Selection  by  FPOD, 
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I’ercent  l’rob;il)ilit  .■  ol'  Detect i(in  lor  Scenario  .'i 
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la  th(  nnifortn  ease  I'esaUs  in  Table  show  that  there  are  still  significant  differences 
in  mode!  mapiiis  at  all  dcptlis.  However,  shoi-tening  the  on  station  time  appears  to  have 
broughi  the  values  closer  together  in  the  normal  case.  'There  is  only  2.3'''  maximum  dif- 
fevence  at  the  S-S  deptii. 


Scenario  I 

1  Te  circle  and  ehev won  skew  iAppendix  Ai  buoy  fields  were  considered.  I  lie  environ- 
meni  was  c  langed  from  the  pacific  to  the  Atl;uitic  for  a  different  type  of  compai’ison.  The 
socii.i  speed  p’rofiles  in  I  igures  li  and  •'!  show  tlia'  there  is  no  CZ  so  sound  propagation  is 
l)V'  (iireel  ivitli.  'I  he  i'roixig;r, ion  loss  curves  in  Ap)ien  lix  IT  indicate  itoltom  Irounce  is  pre.-en: 
wiil'  a  r  -ion  ft  .^oin-ec  rceeiver  <k'pih  comliinat ion.  \'alues  in  I’able  10  sliow  average  Pi  d) 
’  alucs  I'i'taincd  'or  the  Atlantic  ein  i  ronment .  Ilvdrophone  depths  of  100  and  inuO  ft  were 
clio.>e-r  foi  siriulation  and  the  results  ait  shown  in  'I'airlc  10. 

T'l'  -  enaiio  na'  ,  hesen  in  show  l.ow  llie  niodels  respond  in  a  d.ieci  ptilli  eiiv  i  roiimetu , 
A-.  (  co  d,  i(  .all  I  ^  ai'j  I  liver  than  in  the  other  seenanios  because  of  the  tibseuce  of  a 

eot'ver'  '  ni  <■  /ni  c  en  i  r<  ainit  :r  .  the  one  noticeable  result  in  the  output  averaged  in  Table 

1()  i.  pH  t.ijh  ev  dilation  In  .\'/(  i|  using  the  chevron  skew  I'aitern  in  the  normal  case.  .Actu- 

allv  ,  t'l]  .  I  11  '■  I  one  >.'.111111  j a  cdi<  t  for  the  chevron  .Pcevv  (lattern  since  the  coverage  is 

gTt  atcsi  It  I  h(  c'  litcr  n  here  t lie  target  is  most  likely,  i  he  other  two  models,  however, 
did  not  a  lea^t  Pici  1  ■  1  c' iic;  ion  -  -ignificamh  from  I  lie  uniform  case. 


rnl)lc  10 

I'orceni  I’rol  abilitv  of  Oelfction  for  Scenario  4 
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\  !.  CONf  I.CSK  i\S 

I  xpectei!  i.)rol)al)ilit\'  of  detection  is  calculated  using  four  different  methods.  ADId’S 
calculates  the  ttverage  conditional  probability  in  a  symmetrical  subsection  of  a  buoy  field 
ptutern  where  the  sul)marine  is  assumed  to  be  a  random  target,  'I'ASDA  computes  prob- 
abilii\  of  detection  br  one  or  more  buor.s  using  Monte  Carlo  simulation  of  100  random 
tracks  through  a  buov  field  and  a  detect  (POD--  li  or  no  detect  IPOD  -  Oi  analvsis  for 
t  acli  track.  A'/J  d  uses  tt  binominal  probability  function  to  calculate  KPOl)  where  a  \  alue 
lor  lime  along  a  track  is  eho^,en  randomly.  SPAM  uses  a  Monte  Carlo  simulation  and 
c.aleul.ales  )  pi  iJt  bt-  atcr, aging  a  probalnlity  obtained  through  assuming  independence  of 
liuors  ■..■itlt  a  iiro!iabilit\'  oltiained  1)\'  assuming  complete  dependence  of  buoys,  since  the 
m(  iho(U  of  calculation  tire  different  in  each  program,  the  output  for  IfPCD  is  not  exactly 
tlie  same. 

ADI  I’s,  l  AsDA,  A/a  i|  and  SPAM  all  calculate  cumulati\  e  probabilila  of  detection  as 
a  function  o!  time.  However,  llie  method  of  calculation  is  different  in  all  four  programs. 
ADI  PS  u.-^es  an  rxponential  exiiression  to  estimate  a  conditional  cumulative  probability 
I’lti  as  (lisC'u>sod  in  Section  1,  l  ASDA  uses  a  detect  -  no  detect  Monte  Carlo  methodologv. 
A/,(  >1  uses  I  !’(  il)  calculation  to  olittun  an  estimate  for  the  lower  limit  of  CPOD.  In  the 
case  where  the  target  location  is  assunted  to  have  a  uniform  distriltulion,  the  CPCD  does 
not  c  hange  with  time  and  is  equal  to  I  pctD,  In  the  normal  case,  PPOD  and  CPCtl)  are  dif¬ 
ferent  since  ctdculalions  for  the  normal  case  are  time  dependent.  SPAM  uses  a  cumulative 
lime  window  to  limit  l’(  iD  calculations  for  CP()D.  Since  different  techniciues  arc  used  for 
each  program,  it  is  not  eN|iected  that  Cl’ctD  output  will  l)e  exactly  the  same  for  each  pro¬ 
gram.  ADI  I’S  outputs  P  attd  P(t);  1  ASDA  out|)uls  CPCl),  PPCD  and  mean  lime  to  first 
detection  (M  l  i  l>i;  A/(  H  outputs  CI>()D,  I  POD,  sttindard  deviation  (SIC.MA)  and  proltabil- 
itv  tliai  a  target  will  Ite  detected  given  that  it  h;ts  not  Iteen  detected  in  the  past  (!'  Dl  'NDI'i; 
spAM  outputs  Cl'cil)  and  I  pctl).  Although  tlte  svmitols  Cl'OD  and  IfpttD  in  ftible  <>  are 
used  lot  1  ASDA,  A/.(  d,  and  SI'AM,  the  output  for  these  (tarameters  ntav  not  be  the  sante. 


n(  l(  1 1  !i  .'A  ini;  nt'i  iiKlir.'iUfi  I)'.  I  hi "liii.U. 


!.  1  :ili\ (-•  at  I.dw  I't  iMs 

liililsA  .\/(  1]  Of  ''I’AM,  I  ASDA  i'C'<|Ui  fo.-,  coptaoi  (hiring  luo  coiisfi'iitiN  e  time  iteii-i 
l)L  !(irA  a  lit  tfA  tion  is  n.aile,  1  .\Sl)A  i-odst fiu't  s  |iroli:il)ilit\  iii’e.'i.s  afcjuncl  each  l)iio\  Ixmiideil 
1)\  i"i.  foiu'i.  lit  f  ie  t'irc'les  :is  if  |  i^ina  U>.  Iluse  ai'eiM  feprcsent  tlie  loenlioiis  wliei'e  the 


.  IIAAE  STEPS 
i_  i  e  I  1  1  I  1  I  1  1  I  I 

tARGEt  TRACK 


I  |l<51 

CC'NVI  it(.LN(.l 
/ONE 


tASDAIJI  11  1  riON  t'OlNTS 

.  i  t  J 

1  1  j  4  ^  -  * 


DIRECt 


PATH 


i  ieiife  Ih.  Taf'/ei  I'l  aok  fof  l  ASD.X  and  Sl\  >n  ia'i;tor,  Ai'onnd  a  Sonoljuov,  l>(  il)  -  1. 

siunal  excess  SI'  is  jircatcf  than  zero.  If  the  laraei  la  iiiains  in  one  of  those  areas  diii’in'i 
‘.'.'•o  eonstcutix  c  time  steits,  tliontlie  !'(>l)  1  for  ihtit  t  i  me  i>erio(i.  i  ilherwise,  t  lie  !’(>!)  = 

f'.  'I'here  A-t’l  iio  a  set  of  these  \  alnes  foi-  etich  of  the  100  '|  ASDA  tritil.s.  I  he  \  tiliies  are 
sumnied  o'  ei  100  trials  to  ei.inu  op  nitli  Cl'tiD  '.'tiUn.-.s  for  etich  time  interval.  At  low  I'tiMs 
the  aretis  hounded  hv  concentrie  eircUrs  liecomc  vto'v  smtill  aiifl  as  u  result,  the  ('I’f)D 
values  drop  rapidlv  to  zero,  A'/.(  il  and  SDAM  cismputc  I’tiD  versus  ran!i,e  let'  each  Ijttov 
i-atht  r  than  assitininti  ;i  0  or  1  \ nine  tit  a  particular  rajrre. 

'J.  SP.AOI  is  I  n.itimistic  at  Low  I'liMs 

If  the  tarj/et  is  tit  a  raniio  .^uch  that  a  faint  line  will  be  created  on  a  ttram,  then 
tiiei'C  trust  he  some  tin.ount  of  time  durin*;'  which  the  tareet  remains  at  a  detectahle  ratv-tc 
in  order  tor  the  operator  lo  see  a  line:  ctill  it  tf),  !■  \  en  after  tQ  time  has  passed,  the  oper¬ 
ator  mtiv  not  see  the  line.  I'he  pi’ohtihilitv  1  hat  the  operttlor  sees  the  line  inereases  up 
untd  -.I  no,-  time,  tj.,  I  hc'-e  fwii  v.-due.--  fot  time  tito'  incoT'porttfed  into  a  cumultiiiv  e  time 
..'indo'v'  t(  'i  W  I  model  which  ex.nnines  the  proha  hi  I  i  ties  between  these  limes  in  inoT-ements 
of  time  1 1.  I  he  smaller  the  \  alue  of  t(),  the  hijeheT- '  he  l'()I)  value.  Slantlard  nperaline, 
values  !'o;  s|>  \vi  .n-e  ttp,  tl,  'y>  ~  t.-  ',  .AOi.  This  means  that  the  larpct  must  re¬ 

main  at  a  detectahle  ranve  Lo-  onlv  lA  minutt-s  lo  .achieve  the  probnbil  tv  o'  detection  pi'c  - 
dieted  hv  tic  standard  iiassi' c  -.on.a)' efjUation.  The  valnos  tp,  i^,  i;.  in  the  ('TW  eaused 
pt  tdielio  s  ''oi'  SPAM  to  ije  oiitimisiic  at  low  I-'ttAIs.  Detection  can  occur  at  anv  lime 
helorc  or  after  If),  dei-cndino  on  cliance, 

A,  SPAM  ftesults  ;tre  more  Con.scrv ntiv e  at  lli!!:h  l'(  iMs 

I  ne  <  IA\  model  builds  dtnection  probabilities  experimentally  for  a  target  re- 
rnaiiiin  ;  .at  a  fixeii  '  anui  an'd  IItik  '7  vvilli  the  iirohability  pT  t'dicterl  hv  the  sonar  ei|ualion 
atdiie  •  !  I'tej-  titiu'  i,,  has  pa  i.-ed.  )his  funelion  aiso  smoothes  the  lump  from  low  proh- 
ahilitii.  a-si  ciatevl  witii  low  !-(P.i.s  to  hiili  proh.dflilies  asscicialetl  with  hi^h  LtiMs. 


■1 .  I  \SI)A  and  A7,(  )I  Show  a  lump  in  CI’OI)  as  FOM  Inc  reases 


A  normal  distnhution  around  the  fi^aire  of  merit  is  used  in  the  SPAM  model, 
riu'  iielief  is  that  the  area  under  a  normal  curve  is  equal  to  the  probaiDility  of  detection. 

SI  =  0  when  FtiM  -  PI..  Al  that  point  there  is  a  30'','  POD.  As  shown  in  Figure  17, 

SPAM  P(d)  graduallv  increases  as  SF  increases.  Truncation  has  an  effect  on  CPf)I)  values 
for  .\7.<  with  low  Ft'tMs.  For  low  FfiMs  there  is  verv  little  signal  excess  and  prohahilitv 
dei-iendcd  on  the  left  half  of  the  cur  ves  above;  therefore,  the  CPDD  values  dropperl  cjuickh' 
to  I).  IWSDA  uses  a  Pician  Distribution  to  simulate  short-term  fluctuation  in  the  PL  curve 
■  ind  calculatiirns  are  based  on  the  assumption  that  below  the  Ff)M  an  operator  is  unlikely 
to  detect  a  target.  Detection  is  either  VfkS  or  NO  based  on  SF  0  after  employing  the  Pit  iar. 
Huctuations.  As  is  apparent  from  the  graphs  in  Figures  17,  18  and  19.  AZOI  and  TASDA 
ha\  e  probal'ility  distributions  which  drop  off  ciuickly  to  zero  below  the  FONI.  Therefore, 
r.\SDA  and  AZOI  exhibit  breaking  points  when  the  FOM  reaches  a  small  enough  v  alue.  Al 
these  FOMs  (  POD  falls  rapidly  to  zei'O. 


■),  AZOI  Shows  Snrall  I'ariairce  Between  Patterns 

AZOI  differs  from  TASDA  and  SPAM  in  that  TASDA  and  SPAKI  use  successive  time 
Steps  along  a  target  track  for  a  particular  trial.  After  a  track  has  been  generated  for  a 
I)articular  trial,  AZOI  picks  random  points  along  that  track  and  evaluates  POD  values  at 
those  points.  Therefore,  the  time  increments  are  not  evenly  spaced  but  randoml}’  spaced. 

.\t  each  ixDint,  the  range  to  each  buoy  is  calculated  and  the  px’obabilities  associated  with 
each  buov  are  summed  to  give  a  POD  for  the  target  at  a  point  along  the  track.  In  this 
manner,  all  contriixutions  arc  considered.  The  probability  at  a  point  along  the  track  is  not 
zero  or  one  as  with  TASDA,  nor  is  execution  terminated  when  a  detection  is  encountered  as 
with  SPAM.  So  this  method  yields  more  values  for  POD  along  a  track;  and  a  cumulative 
time  window  is  not  employed  as  is  the  case  with  TASDA  and  SPAM. 

llicse  five  differences  in  the  three  computer  simulation  models  result  from  various 
interpretations  of  the  tactical  sonobuoy  situation  under  consideration.  These  interpretations 
result  in  (.lifferent  output  values  for  each  of  the  three  simulation  programs. 

I'he  large  differences  between  computer  models  require  further  investigation.  Future 
Studies  should  ije  based  on  a  real  world  application  of  the  sonobuoy  pi’ograms.  Ileal  world 
mission  planner.s  usuallv  have  SOSFS  or  other  intelligence  information  for  input  to  the 
model  and  are  intei  ested  in  POD  versus  buoy  spacing.  The  output  data  sets  in  this  report 
are  not  designed  to  i'eflect  a  real  world  application  of  the  sonobuoy  pi’Ograms,  but  acquaint 
the  reader  with  the-  amount  and  type  of  data  output  generated  by  considering  a  few  pattern 
cfintigu I'ations  under  several  different  conditions. 
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Figure  A-1.  Circle  pattern  used  in  scenarios  1,  2,  3,  and  4 
for  a  uniform  target  distribution. 
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Figure  A'2,  Ellipse  pattern  used  in  scenarios  1,  2,  3,  :md  4 
for  a  normal  target  distribution. 

A-" 


Figure  A-3.  4-4-4-4  pattern  used  in  scenario  1  for  a  uniform  and 
normal  target  location  distribution  and  in  scenarios  2, 
3  and  4  for  a  uniform  distribution. 
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Figure  A-4.  Chevron  or  X  pattern  used  in  scenario  1 
for  a  uniform  target  distribution. 


i  inin'o  A-.').  Chevron  skew  pattern  used  in  scenarios  1,  2,  2,  and  4 
tor  a  nornud  target  distribution. 


{ _ 

l■■:g^Jre  A-G.  4-0-5  pattern  used  in  scenario  1  tor  a  uniiorm 
and  normal  target  location  tlistribuUon. 
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riguro  A-7.  Brushtac  pattern  used  in  scenario  1  for  a  uni¬ 
form  target  distribution  and  for  a  normal 
target  distribution. 


A- 5 


A!>[>i;ni)ix  li 


PUOI'AGATIOX  LOSS  CIKVLS 


B-1 


tsl  H  H 


-J  O  tf) 


‘7:^  S  10  15  £<?  25  30  35  4©  4S  5S  6«  65  7®  7S  80  8S  90  96 

KHf-tCitTl _ _  .. — . — —  -  —  - 


RANGE  KVDS 

FREQUENCV . 

1000  HZ 

80  1 

00  i 

80  140  160 

180  20 

■■ 

■■ 

■ 

■ 

B 

B 

Hi  Hi 

mnn 

n 

IB 

Bl 

■ 

Bl 

B 

B 

IB 

HIM 

HHI 

IH 

■ 

Bi 

B 

B 

B 

MHi 

■■ 

n 

H 

■ 

B 

B 

B 

B 

■n 

»^Hi 

[B 

■ 

B 

B 

B 

B 

22 

SUIS 

9 

2 

2 

51 

m 

2 

iP 

■i 

2 

2 

■■1 

HIHI 

Hi 

Hi 

■ 

■ 

s 

91 

B 

i5 

5  10  15  30  as  30  35  40  45  50  55  60  65  70  75  TO  TO  90  96 


I'ifliire  B-2.  Propagation  lows  curves  at  a  sh;dlo\v  receiver  depth  mm  deep 
source  depth  for  scenarios  1,  2,  and  3. 


n.  ii:  <:>  ri,  _i 


PPOPAGAT I ON-LOSS  GRAPHICS 

DATE  -  7  '  6  / 

LAT  -  3558  N 

LONG  -  13559  U 

71 

UAOE  HEIGHT . 

5 

FT 

SOURCE  DEPTH. . . 

30C 

FT 

SONIC  LAVER  DEPTH.. 

FT 

RECEII.CR  depth. 

300 

FT 

FREQUENCV . 

50 

HZ 

PPOf>^«G^TION-LOSS  GRAPHICS 


Figure 


Um'.'E  HEI'GHT .  5  FT 

CONK  LAVER  DEPTH..  lia  FT 


DATF  -  11  18  /■  69 

LAT  -  393S  M 
LONG  •  7130  u 


SOURCE  DCPTH. . .  15«  FT 
RECEIOER  depth.  100  FT 
FREOUENCV .  100  HZ 


PPOPAGATION-LOSS  GRAPHICS  'FNUC  9  ?•»  ' 


UAiJE  HEIGHT .  5  FT 

SONIC  Laver  depth..  iie  ft 


DATE  •11/18 
LAT  .  3935  N 
LONG  •  7130  U 

SOURCE  depth.  .  . 
RECEIVER  DEPTH. 
FREQUENCY . 


69 


500  FT 
1000  FT 
100  HZ 


B-4.  I’ropagation  lo.ss  curves  exhibiting  direct  path  sound  tnmsmission 
with  no  convergence  zone  used  as  input  to  the  tactical  models  in 
scenario  4. 
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C  43.  Comparison  of  TASDA.  A/OI  and  SPAM  witli  a  uniform  initial 
lai’irf't  location  at  one  liou>’  tinm  late  with  a  shallow  taru'ct 

.i  hallow  receiver  i.i  .axai  irif;  ,'i  n.siiig'  the  circle  pattern. 
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SCENARIO  4  -  UNIFORM  CASE 


t  i  .i.iri'i'i,  .)|  lASDA,  AZOI  .  .iiut  SPAM  witli  a  uniroriii  initial 
U>ka;i''::  at  o.’u-  lu>i:r  ti:iio  lato  with  a  ix'oeivor  li'i’quoncy 
'■!  aiO  II/.  m  the  (liilf  Stream  on virnnmoiit  in  sconario  1  usinn' 
tiio  eiri  lo  Dattorn. 
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Figure  C-57.  Comparison  of  TASDA,  AZOI  and  SPAM,  with  a  uniform  initial 
larg'd  location  at  one  hour  time  late  with  a  receiver  frociucncy 
of  100  Hz  in  the  slope  water  environment  in  scenario  4  using 
the  circle  pattern. 
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Figure  C  58.  Compni-ison  of  f  ASDA .  .\7.0\  and  SPAM,  with  a  viiforiii  initial 
tariTot  location  at  one  hour  time  late  with  a  receiver  fre<)ucncy 
of  >100  H/  in  the  slope  water  environment  in  scenario  i  using 
the  circle  pattern. 
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Figure  C-59.  Comparison  of  TASDA,  AZOI ,  ami  SPAM  with  a  itoriiial  initial 
target  location  at  one  hour  time  late  witli  a  roci'ivcr  l'ri'i|ucnc\ 
of  100  Hz  in  the  (iulf  Stream  environment  in  scenario  -1  using 
the  chevron  pattern. 
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Figure  C-GO.  Copip.irisoii  of  TASDA ,  AZOI.  :in  SPAM  witli 
tr.i’LV'  t  locTitioii  at  ono  hour  tiiiu'  latf  u'it)i  a  r. 
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P'igure  C-03.  Scenario  1  -  Fventual  I’robability  of  Detection  for  a 

shallow  target  and  receiver.  The  eventual  I'OD  value 
is  tl  '■  maximum  ''i  I’OD  attained  by  the  cur\'e.  i:ventual 
POD  is  obtained  by  obser\ing  the  maximum  POD  for 
each  CPOD  curve.  These  maximum  values  are  summed 
over  frequency  and  FOM  to  give  a  percent  POD  for  a 
uniformly  distributee!  shallow  target  and  a  shallow 
receiver.  Not  for  example;  (.  14  +  .50  +  .90  j  .22  + 

.58  +  .86  *  .21  +  .26  ^  .76)/9  =  49.9. 
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M  l'ri)l)nl)ilii\  of  Detoi  l ion 


Table  C'-3:  Scenario  I 
Kventual  Probability  of  Detection 
for  i  I'niform  'rart*ct 

at  30  min  Time  Tate:  Shallow  -  Deep  Case 
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Table  C-9:  Scenario  1 
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I’ablc  C'-12:  Scenario  1 
Eventual  Probability  of  Detection 
for  a  Normal  Target 
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Tal3let''13:  scennrio  2 
Average  Probability  of  Detection 
for  a  Uniform  Target 
at  1  Hr  Time  Late 
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l  abel  (■- 11:  Seenario  2 
Average  l*n)bability  of  Detection 
for  a  Normal  I’arset 
at  1  Hr  Time  Late 
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Figure  C-65.  Scenario  3  and  4:  Average  Probability  of  Detection  for  a 
shallow  target  and  receiver.  The  average  POD  value  is 
obtained  by  summing  the  POD  values  at  each  time  incre¬ 
ment  .and  dividing  by  the  total  number  of  time  increments 
to  estimate  the  area  under  the  CPOD  curve.  The  average 
values  are  summed  over  frequency  and  FOM  to  give  a 
percent  POD  for  a  uniformly  distributed  shallow  target 
and  a  shallow  receiver. 
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111.0 

1520.5 

13a. 0 

1519.2 

147.0 

1513.5 

165.0 

1516.7 

181.0 

1515.8 

200.0 

1514.5 

312.0 

1501.7 

460.0 

1491.2 

600.0 

1485.5 

800.0 

1481.7 

1000.0 

1481.4 

1200.0 

1483.8 

1500.0 

1488.9 

2000.0 

1497.4 

2400.0 

1503.5 

SOUND  UELOCITY  PROFILE  LAT  3936N  LON  7029U 

UELOCITY  n/SEC  DATE  171169 


